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Abstract 

RICH detectors have become extraordinarily useful. Results include measurement of 
solar neutrino rates, evidence for neutrino oscillations, measurement of TeV 7-rays 
from gravitational sources, properties of QCD, charm production and decay, and 
measurement of the CKM matrix elements V cs , V c b and V u b- A new value |T4b/V c fe| = 
0.087 ± 0.012 is determined. 
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1 Introduction 

This paper describes some physics results obtained by experiments using Ring 
Imaging Cherenkov Detectors. Currently, the major physics issues are: 

• What is the origin of mass? Another way of phrasing this question is to state 
that the Higgs boson must be found. This is a crucial reason for building the 
LHC machine. In this case, the RICH detectors have not yet been shown to 
have any relevance. 

• Do neutrinos have mass? If so can neutrino mixing be explained in the 
context of the Standard Model? 

• Does Quantum Chromo-Dynamics explain the strong interactions? 

• How to unify gravity with the other interactions? 

• Does the Standard Model explain quark mixing via the CKM matrix? We 
need to measure the CKM elements and CP violating angles. 

The raison d'etre for RICH detectors is to answer these and other questions. 
The measurements reported here reflect my own view and may be incomplete. 

RICH detectors can be viewed as being used for two distinct functions. One 
is to detect the presence of charged particles and the other is to identify the 
kinds of particles that have been detected by other devices. 



2 QCD Results 

One of the first RICH detectors was used in Fermilab experiment E605. They 
measured particle spectra produced in 800 GeV/c proton interactions with Be 
and W nuclei. The detector used a He gas radiator and a He-TEA photon de- 
tector. On average only 3 photo-electrons per track were detected. Yet results 
were produced [1]. 

More sophisticated modern versions of fixed target spectrometers are the 
CERN OMEGA spectrometer and SELEX at Fermilab [3]. A diagram of SE- 
LEX is shown in Fig. 1. Precision silicon strip detectors are interspersed with 
the target. A gas RICH detector with phototube readout is used to indentify 
charged hadrons. 

SELEX will study the decay of charmed baryons. They have already produced 
results on the production ratios of charmed mesons and baryons from different 
beams [4], shown in Fig. 2. These can be compared to theoretical models of 
charm production. 

Two machines, LEP and SLC produce Z° bosons in e + e" collisions. Physics 
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Fig. 1. Diagram of the Selex detector. The dashed arrow points to the vertex and 
target region, that contains the silicon strip modules. 
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Fig. 2. Invariant mass plots for A+ and D + and their antiparticles for different 
incidents beams. The measured production ratios are shown on the figure. 

results have been obtained by two experiments using quite similar RICH de- 
tectors, DELPHI [5] at LEP and SLD [6] at SLC. A primitive sketch of the 
RICH systems is shown in Fig. 3. These systems use both liquid C 6 F 14 and 
gaseous fluorine radiators. The Cherenkov photons are converted in TMAE 
and drifted using a TPC to proportional wires. The performance of the SLD 
CRID has been characterized in terms of efficiency and rejection for a par- 
ticular set of analysis criterion. Shown in Fig. 4 (left) are efficiencies for the 
wanted hadron and the efficiency for the unwanted species (really rejection) .0 

The Z° decays into quark- ant iquark pairs or lepton-antilepton pairs. In the 
case of the quarks, the energy ends up in hadrons and the momentum spectrum 



2 In general, the efficiency versus rejection is a function of the analysis cuts and 
can be chosen differently to optimize each study. 
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Fig. 4. (Left) Efficiencies for tracks from K° (squares) and r (circles) decays to 
be identified as each hadron species in the SLAC CRID. The solid symbols show 
the data and the open symbols the simulation. (Right) Comparison of measured 
charged hadron fractions (symbols) from SLD compared with the predictions of 
various fragmentation models. 

of the hadron species can be compared with models of this fragmentation. The 
particle spectrum from Z° decays has been measured by both DELPHI and 
SLD. The SLD spectrum is shown in Fig. 4 (right) and compares well with the 
different models [7]. 



QCD predicts the presence of a new state of matter, a so called "quark-gluon 
plasma," that may be produced in the collisions of heavy nuclei. Evidence 
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could come dileptons produced in elementary collisions such as qq — > 
Electron pairs are typically searched for. The can also come from more mun- 
dane sources, such as vector meson decays or Dalitz decays. 

The CERES experiment has searched for e + e~ pairs in heavy ion collisions 
at the CERN SPS. A sketch of the detector is shown in Fig. 5(left). They 
have two RICH detectors using a CH4 gas radiator and UV (TMAE based) 
photon detectors, which are placed outside the region of the collision products 
[8]. Also, heavier particles, such as pions, tend not to radiate, so this type of 
system is sometimes called "hadron blind." 




m ee (GeV/c 2 ) 



Fig. 5. (Left) Sketch of the CERES detector. (Right) The cross-section for e + e~ 
pairs as a function of dilepton mass. The yield from expected sources is also shown. 
The statistical errors (bars) and the systematic errors (brackets) are plotted inde- 
pendently of each other. 



Fig. 5(right) shows there results for Pb-Au collisions [9]. There is an excess of 
e + e~ pairs above the estimated background sources. Many new experiments 
are being built to explore this phenomena including HADES at GSI Darm- 
stadt, and RHIC experiments at Brookhaven. 



3 Water Cherenkov Detectors 



3. 1 Introduction 



The quest for a theory to unify strong and electroweak interactions led to 
models based on higher symmetry groups which predicted or at least allowed 
the possibility of proton decay. The 1MB and Kamiokande detectors were large 
tanks of pure water placed deep underground to minimize the cosmic ray 
background. The water provided the protons that might decay and also the 
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means of detection using the Cherenkov light produced by the decay products. 
Kamiokande was greatly enlarged to Super-Kamiokande [12]. 

All these detectors have many large diameter photomultiplier tubes which 
gather the unfocused light (proximity focused in Tom Ypsilantis terminology). 
They can detect charged leptons and hadrons and also photons that convert 
to pairs in the water. They can discriminate between electrons and muons 
because the electrons scatter more. Fig 6 shows the phototubes hit by an 
electron and muon. The width of the Cherenkov ring is much smaller in the 
muon case. 



• Single Electron: • Singfe Muon: 




Fig. 6. Examples of the detector response of an electron and a muon in Su- 
per-Kamiokande. 

These detectors can also detect v e from the sun and from supernovas. The 
latter capability was impressively demonstrated by detecting the explosion of 
supernova 1987A [13]. They also can detect neutrinos produced in the atmo- 
sphere from cosmic ray interactions. 

3.2 The Atmospheric Neutrino Anomaly 

Cosmic rays, mostly protons, hit the atmosphere and interact at altitudes 
of ten to twenty kilometers above the Earth's surface. Naively, we expect 
two muon- neutrinos, v^, for every v e , because of the dominant decay chains 
ir + (K + ) — > /i + — > e + u e i7^ (and similarly for 7r~(K~)). 1MB first observed 

that the ratio was about half of what was expected [14]. The data from several 
experiments are shown in Fig. 7. 

Although some experiments did not see the anomaly, the precise data from 
Kamiokande and Super-Kamiokande [16] clearly show that 1MB was correct. 
The question remained, however, was there something wrong in the expecta- 
tion or was the effect coming from a loss of due to neutrino oscillations? 

In neutrino oscillations, different species can metamorphose into other ones. 
We have three known species of neutrinos, u e , v„ and v T . If u' s have mass 
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Fig. 7. Summary of measurements of the atmospheric neutrino ratio of measured 
muon/electron neutrino rates divided by the expectation as given by Monte Carlo 
prediction. (v^/i'e)data/(i / [_i/ve)MC for atmospheric experiments. The three newest 
results are listed at the right (from ref. [15]). 

they can transform from one to another. Specifically, the probability that a 
neutrino, z/j, not oscillate is given by 

Pfa -». Vi ) = 1 - sin 2 (2fl) sin 2 (l.27Am 2 J^ , (1) 



where 9 is the unknown mixing angle, Am is the unknown mass difference, L 
is the distance from the production point and E v the energy. 

Cosmic rays hitting the upper atmosphere are distributed uniformly. The ex- 
periment is at a fixed point close to the Earth's surface. If a neutrino comes 
from directly overhead the cosine of its angle is denned as +1, while those com- 
ing through the earth have a cosine of -1. Thus u's have different path lengths 
to get to the detector; since the v energy spectrum is the same from any di- 
rection, we can vary the ratio L/E v by looking at different angles. These u's 
range from below 1 GeV to many GeV in energy. (Super-Kamiokande defines 
sub-GeV to be <1.3 GeV, and multi-GeV to be larger.) 

Super-Kamiokande is a large enough detector to get enough atmospheric neu- 
trino events to be able to show angular distributions. While u e behave as 
expected, the distribution of u^ do not. In Fig. 8 (left) the zenith angle distri- 
bution of Up is shown along with the expectations based on simulation. There 
is a clear deficit of up-going neutrinos. On the right side the L/E v distribution 
is plotted for both u e and u^. The evidence for u^ oscillating into something 
is clear for u^. The dashed lines show the expected shape for u^ u T us- 
ing Am 2 = 2.2 x 10~ 3 eV 2 and sin 2 (2#) = l. The dashed lines for u e show the 
expected shape without oscillations. The slight L/E v dependence for e-like 
events is due to contamination (2-7%) of u^. 

If the are oscillating, they can change into u e , v T or some other species of 



7 



yet unobserved neutrino. As the latter possibility is less likely, let us ignore 
it here. We do not see an increase in v e , so its most likely that the are 
changing into v T . Unfortunately, the v T are not energetic enough to produce 
t~ via a charged current interaction, most of the time. 
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Fig. 8. (Left) Zenith angle distributions of multi-GeV events. Upward-going par- 
ticles have cos# < and downward-going particles have cos# > 0. The shaded 
region shows the Monte Carlo expectation for no oscillations normalized to the 
data live-time with statistical errors. (Right) The ratio of data to Monte Carlo 
versus reconstructed L/E v (fully contained events only) [16]. 



3.3 The Solar Neutrino Anomaly 



Thermonuclear reactions in the center of the sun are a prolific source of v e . 
Detecting them and comparing the rate with models describing the process 
have gone on since 1964, starting with Ray Davis' pioneering experiment in 
the Homestake mine [17]. Confirmation that the rate is about half what is 
expected was made by the Gallium based experiments (Gallex and Sage) [18] 
and the water Cherenkov detectors [19]. The expected flux, and the results 
are shown in Fig. 9. 

Super-Kamiokande can also show that the detected v e come from the sun. 
Fig 10 shows the angle of the produced electron with respect to the solar 
direction [19]. The sharp peak is evident above the flat background. 

One explanation of the v e deficit is that the sun is just at the right distance 
from the earth that half have oscillated to other species. This is called the "Just 
So," solution [20]. Another explanation is there is a resonant enhancement of 
oscillations due to the high electron density in the suns core; this is called the 
MSW effect [21]. 
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Fig. 9. (Left) The predicted solar neutrino energy spectrum. The nuclear reactions 
are indicated. (Right) The measured solar neutrino flux from the different experi- 
ments divided by the expectations of the standard solar flux model. Gallium exper- 
iments are sensitive above 0.23 MeV, CI above 0.8 MeV, and H2O above 6.5 MeV, 
from [15]. 
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Fig. 10. The solar neutrino spectrum from Super-Kamiokande as a function of the 
angle with respect to the sun. 



3.4 Summary of Neutrino Oscillation Signals and Future Prospects 



Possible signals for neutrino oscillations come from three sources, two of which, 
atmospheric neutrinos and solar neutrinos have been described here. A third 
source, the LSND experiment claims observation of a signal for — > u e [22]. 
Putting these on the same Am 2 — sin 2 (26 l ) plot (Fig. 11) shows that it will be 
difficult to reconcile all three measurements. 

In the near term future there will be several new experiments using Water 
Cherenkov detectors that should reveal much about the mysteries just begin- 
ning to be revealed. A beam of will be directed at Super-Kamiokande [23]. 
The SNO detector, a large D2O detector will be able to detect the neutral cur- 
rent neutrino interactions and thus be able to measure the total atmospheric 
neutrino flux and the total solar flux as well as charged current interactions 
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Fig. 11. Allowed regions for three indications of neutrino oscillations. The top blob 
for Atmospheric is from Kamiokande while the bottom is for Super-Kamiokande 
(from ref. [15]). 

[24]. Borexino will start and be able to detect 0.86 MeV neutrinos from solar 
Be 7 reactions at the level of 50 events/day [25]. 



4 Measurement of Primary Cosmic Rays with Caprice 

Analysis of the atmospheric neutrino data relies on knowledge of the primary 
cosmic ray flux. The Caprice experiment has taken data with a magnetic 
spectrometer incorporating a RICH detector with a NaF radiator and photon 
detector using a wire-proportional chamber filled with TMAE and ethane [26] . 
They also have tracking and an electromagnetic calorimeter. 

CAPRICE has measured the rates of many species at the top of the atmo- 
sphere. The response to positive tracks is shown in Fig. 12(left). Their primary 
cosmic ray (mostly proton) flux measurement [27] is compared with others in 
Fig. 12 (right). The data that shows a larger rate is from older measurements. 
They have also measured the p/p rate to be ~ 10~ 4 , /i + //i~=1.64±0.08 and 
e + /e~ to be ~ 0.1 

Calculations of the primary cosmic ray flux and their interactions with air, 
producing secondary mesons that then decay into leptons and neutrinos are 
important in evaluating the total neutrino yield and less consequential in eval- 
uating the angular distribution of the z/ M /z/ e ratio. Super-Kamiokande uses 
several calculations of primary cosmic ray flux [29] [30] . The HKKM [29] cal- 
culation of the primary cosmic ray flux is shown as "Mid" on Fig. 12 (right). 
Francke has argued at this meeting that the HKKM calculation over estimates 
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^ Energy per nucleon (GeV) 

Fig. 12. (Left) Cherenkov angle as a function of rigidity, R, for positive particles in 
CAPRICE, where R = (A/Q)(E 2 + 2M Q E) 1 / 2 , M o =0.9315 GeV, A is atomic mass 
number, and Q is charge. The bands, top to bottom, correspond to p and He 
nuclei. (Right) Measured primary cosmic ray proton flux. The CAPRICE data are 
shown by open vertical diamonds [27]. Other data are described in ref. [28]. The 
labels High, Mid and Low refer to the ranges of fluxes used by HKKM [28]. 

the CAPRICE data by about a factor of two. This could lead to some reinter- 
pretation of the atmospheric neutrino data and may change the limits in the 
Am 2 — sin 2 (2$) plane. 



5 Gravitational Physics with Cherenkov Detectors 

The Fly's Eye detector in Utah uses scintillation light to detect ultra-relativistic 
cosmics rays [31]. Cherenkov light produced in the atmosphere from TeV 7- 
rays, after converting to e + e~ pairs can also be used. HEGRA incorporates 
several Cherenkov light detectors with scintillation detection of electromag- 
netic particles, and muon detection over an area of 200x200 m 2 [32] Q 

One example of what they have observed is shown in Fig. 13. They measured a 
high flux of TeV 7 rays from the source Mkn501. The energy spectrum follows 
a power law from 1-10 TeV. Furthermore, the source intensity varies from 
half that of the Crab Nebula to six times in very small time intervals. This 
phenomena has not been explained, but clearly these measurements will lead 
to an increased understanding of gravitational physics phenomenology [33]. 



6 Measurment Of The Cabibbo-Kobayshi-Maskawa Matrix 



3 For a descripton of other Cherenkov air shower experiments, see E. Lorentz in 
these proceedings. 
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Fig. 13. Detection rate of Mkn 501 on a night by night basis (a) for the whole 
data sert and in 5 min. intervals (b) for the last 3 nights. The dashed lines give the 
average per night, the dotted line gives the Crab rate. 



6. 1 Introduction 



In the "Standard Model" of Electroweak interactions, the gauge bosons, W ± , 
7 and Z° couple to mixtures of the physical d, s and b quark states. This 
mixing is described by the Cabibbo-Kobayashi-Maskawa (CKM) matrix [34]. 



^ Vud Vus Vub ^ 



V ( 



CKM 



Vd V cs V c b 
Vtd V ts Vtb 



(2) 



where the subscripts refer to the quarks. In the Wolfenstein approximation^ 
the matrix is written in terms of the parameters A, A, p and r\ as [35] 



VcKM — 



1 - A 2 /2 A AX 3 (p - irj) 

-A 1 - A 2 /2 AX 2 

AX 3 {1 -p-irj) -AX 2 1 



\ 



(3) 



These parameters are fundamental constants of nature that need to be 
determined experimentally, as is required for other fundamental constant such 
as a or G. 



4 In higher order other terms have an imaginary part; in particular the V c d term 
becomes —A — ^4 2 A 5 (p + irj), which is important for CP violation in decay. 
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Fig. 14. The regions in p — rj space (shaded) consistent with measurements of CP 
violation in iif£ decay (e), V u b/V c b in semileptonic B decay, B° d mixing, and the 
excluded region from limits on B° mixing. The allowed region is defined by the 
overlap of the 3 permitted areas, and is where the apex of the CKM triangle sits. 
The bands represent ±ler errors. The error on the Bd mixing band is dominated 
by the parameter Here the range is taken as 240 > Jb > 160 MeV. |Ktb/K*| is 
taken as 0.087±0.012 (see section 6.3). 



The constants A and A are determined from charged- current weak decays, and 
are 0.22 and ~0.8, respectively [36]. The r] term gives rise to a complex phase 
in the matrix which allows CP violation in weak interactions. A similar matrix 
may exist for u's, if they have mass. It would explain v mixing and predict 
the possibility of CP violation in v interactions. 

Although p and rj have not been determined, there are measurements that 
provide constraints on their value. These include CP violation in K1 decay, 
characterized by the value e, B° ^ B mixing and Vub/Vcb. My interpretation 
of the current status of these measurements is shown in Fig. 14. Here I have 
used only ±ler errors, with the caveat that the dominant errors in all three 
cases are caused by estimates of uncertainties on theoretical parameters. 

Study of these processes has suffered from lack of particle identification ca- 
pabilities. CLEO and ARGUS had time-of- flight scintillators and dE/dx mea- 
surements from the tracking chamber [37]. However, they are blind to K/tt 
separation between ~ 1 — 2.2 GeV/c with poor separation above 2.2 GeV/c. 
Yet CLEO was the first to measure \ V u b/V c b\ and both made the first relatively 
precise measurements of \V c b\. Recently, the DELPHI RICH has brought some 
improvements, which point to future possibilities. 
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6.2 Measurement OfV c b Using B — > D*£is 



Currently, the most favored technique is to measure the decay rate of B — > 
D*£~u at the kinematic point where the D* + is at rest in the B rest frame 
(this is often referred to as maximum q 2 or u = 1). Here, the theoretical 
uncertainties are at a minimum.^] There have been previous results using this 
technique using the decay sequence D* + — > tt + D°; D° — ► K~n + , or similar 
decays of the D*°. 

In a recent analysis, DELPHI detects only the slow 7r + from the D* + decay 
and does not explicitly reconstruct the D° decay [39] . The RICH helps in mea- 
suring missing energy in this complicated measurement. The distribution of 
events with respect to u is shown in Fig. 15. The dip in the data near u of one 
is an artifact caused by the drop in efficiency in detecting the slow 7r + at low 
momentum. Still, the decay rate can be ascertained. Table 1 summaries deter- 
minations of V c b- Here, the first error on is statistical, the second systematic 
and the third, an estimate of the theoretical accuracy in predicting the form- 
factor at uj of 1. Currently, DELPHI has the smallest error, however, CLEO 
has only used 1/6 of their current data and will certainly improve on this mea- 
surement. The quoted average |V^,| = 0.0381 ± 0.0021 combines the averaged 
statistical and systematic errors with the theoretical error in quadrature and 
takes into account the common systematic errors, such as the D* branching 
ratios. 
Table 1 

Modern Determinations of using B — > D*l~V decays at u = 1 

Experiment V c b (xlCT 3 ) 

DELPHI [39] 41. 2 ±1.5 ±1.8 ±1.4 

ALEPH [40] 34.4 ± 1.6 ±2.3 ± 1.4 

OPAL [41] 36.0 ±2.1 ±2.1 ±1.2 

CLEO [42] 39.4 ±2.1 ±2.0 ±1.4 

Average 38.1 ±2.1 



6.3 Measurement OfV u b 

Another important CKM element that can be measured using semileptonic 
decays is V u b- The first measurement of V u b done by CLEO and subsequently 
confirmed by ARGUS, used only leptons which were more energetic than those 

5 Current estimates of the form-factor necessary to translate the decay rate mea- 
surement into a value are 0.91±0.03 [38]. 
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Fig. 15. Fit to the uj distribution from DELPHI. The solid points are the data, the 
light shaded area corresponds to combinatorial background, the dark shaded area 
correspond to other background and the horizontal line contains these components 
plus the signal. 

that could come from b — > ct~v decays [43]. These "endpoint leptons" can 
occur, b — > c background free, at the T(4S) because the S's are almost at 
rest. Unfortunately, there is only a small fraction of the b — > ui~v lepton 
spectrum that can be seen this way, leading to model dependent errors in 
extracting final values. 

DELPHI tries to isolate a class of events where the hadron system associated 
with the lepton is enriched in b — > u and thus depleted in b — > c [44]. They 
define a likelihood that hadron tracks come from b decay by using a large 
number of variables including, vertex information, transverse momentum, not 
being a kaon. Then they require the hadronic mass to be less than 1.6 GeV, 
which greatly reduces b — > c, since a completely reconstructed b — > c decay 
has a mass greater than that of the D (1.83 GeV). They then examine the 
lepton energy distribution for this set of events, shown in Fig. 16. 

DELPHI finds \V ub /V cb \ = 0.104 ± 0.12 ± 0.015 ± 0.009. The first error is 
statistical, the second systematic and the third the uncertainty quoted by 
Uraltsev on his model that allows the extraction of \V ub \ from the measured 
branching ratio [45]. The systematic error is larger than the statistical error. 
This reflects the extensive modeling of the b — > c decays. The theoretical error 
is estimated at 8%. However, another calculation using the same type of model 
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Fig. 16. The lepton energy distribution in the B rest frame from DELPHI. The data 
have been enriched in b —* u events, and the mass of the recoiling hadronic system 
is required to be below 1.6 GeV. The points indicate data, the light shaded region, 
the fitted background and the dark shaded region, the fitted b — > ulv signal. The 
theoretical models are described in ref. [47]. 



by Jin [46] gives a 14% lower value, with a quoted error of 10%. 

ALEPH [48] and L3 [49] have used techniques similar to DELPHI. I have 
averaged all three LEP results and shown them in Fig. 17 without any the- 
oretical error. My best estimate of [K^/V^l using this technique includes a 
14% theoretical error added in quadrature with a common systematic error 
of 14%, since the Monte Carlo calculations at LEP are known to be strongly 
correlated. 

Also shown in Fig. 17 are results from CLEO using the measured the decay 
rates for the exclusive final states niv and p£u [51], and results from endpoint 
leptons, dominated by CLEO II [50]. From the exclusive results, the model of 
Korner and Schuler (KS) is ruled out by the measured ratio of p/n. This model 
deviated the most from the others used to get values of \V u b\ from endpoint 
leptons. Thus the main use of the exclusive final states has been to restrict 
the models. The endpoint lepton results are statistically the most precise. 
Assigning a model dependent error is quite difficult. I somewhat arbitrarily 
have assigned a 14% irreducible systematic error to these models and used 
the average among them to derive a value. My best overall estimate is that 
\V ub /V cb \ = 0.087 ±0.012. 
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Fig. 17. Measurements of using different techniques and theoretical models 

[47]. (The KS model has been ruled out.) 

6.4 Measurement of V cs 



Using the reaction e + e~ — ► W + W" at LEP II DELPHI has directly measured 
the CMK element V cs There are two ways that LEP II experiments have 
estimated this value. One way, is to measure the branching ratio of the W's 
into hadrons and use the sum rule: 



3r^ Br(W ± — > hadrons) 
\\ T h (1 - fir (W ± -»• hadrons)) 



(4) 



where Ti and Th, are the leptonic and hadronic widths of the W 



± 
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A better way, which does not rely on other measurements to measure the 
ratio directly of Br(W — > cs)/Br(W — > hadrons). For this DELPHI uses the 
RICH to identify kaons, which come not only from the s quark but are also 
prolifically produced in c decays. Using this method DELPHI finds 

|vy = 1.01+^ ±0.10, (5) 



which is by far the best direct measurement of this quantity, and is consis- 
tent with the expected value of 0.95. This is important because it checks the 
predictions of the CKM hypothesis.^ 



7 Conclusions 



It takes a Village to build a RICH, a village of knowledge cultivated at this 
conference. There are several different technologies: gas, H 2 0, and crystal ra- 
diators coupled with gas, Csl, phototube photo-electron detectors, producing 
great results in neutrino physics, b decays and QCD. We expect even more 
interesting results soon! 
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6 The indirect method gives (V^l = 0.98 ± 0.07 ± 0.04 
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